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1. Introduction. 

1.1. A Recollection of Preliminary Notions. The Hamilton- Jacobi formula- 
tion of Classical Dynamics is usually presented [2, 3, 10, 11, 16, 18, 26, 32, 41, 47, 
51, 55] as the search, for a Hamiltonian system with a possibly time-dependent 
Hamiltonian on the cotangent bundle T*Q of some configuration manifold Q (with 
dim Q = n for some n) , for a canonical transformation that is able to " reduce the 
system to equilibrium" .The Hamilton- Jacobi (H J) equation for the generator S of 
the transformation, also known as " Hamilton's principal function" , is well known 
to be then: 

dS \ dS 



where: H = H (q;p;t) is the Hamiltonian of the system and, if a complete integral 
S = S (q; Q), i.e. a solution depending on as many additional parameters Q , Q n 
as the number of degrees of freedom in an essential way, i.e. such that: 

det „ ^0 (2) 



dtfdQi 

is available, then, the canonical transformation: (q,p) — > (Q,P) defined by: 

dS dS 
Pi= W , Pi = ~QQV » = !>•••» ( 3 ) 

does the job of reducing the system to equilibrium. 

For a time-independent Hamiltonian, and hence for a conservative system, de- 
noting by E the total energy, the assumption: 

S = W - Et (4) 

yields instead the time-independent HJ equation: 

bUT\-B (5) 



dq 

for "Hamilton's characteristic function" W, with the energy entering as one of the 
additional parameters on which the principal function has to depend. However, 
even in this case, the most general solution of Eq.(l) need not be linear in t. 

The HJ equations come from the search for a canonical transformation on the 
extended phase space R 2n x R such that: 

( Pl dq l - Hdt) - (PidQ 1 - Kdt) = dS (q, Q; t) (6) 

with the additional requirement that K = 0, i.e. that, in the sought-for new 
coordinates: 

or, more generally: 

be the "model dynamics" we would like to relate with our starting one, i.e.: 

d dH dq 1 dH 

Jt Vi ~~~dq~^ ~dt~dp~ 
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Of course, we might consider other "model dynamics". For example, we might 
require that: K = K (P) = (Pj 2 + ... + P^) /2, in such a way that the "model 
dynamics" would be: 

_ (10) 

Clearly, other models could be implemented. For instance, we might require: 

^ / ■ \ 2 2 

and hence: K = *Ylii v A\Q % ) + (Pi) The associated partial differential equa- 
tions would be: 

dS \ dS 1 /dS s 2 



or: 



H ^^^ + m=2(dQ) (12) 



as \ ds v 



(Q J 



( dS 



(13) 



While the first two cases are strictly local, the third one is "less local" , it identifies 
some equilibrium points, say Q l = 0, Pi = 0, which are also stable ones. 

The conditions for solvability of these PDE's are quite strong. For instance, 
the first one requires: dPi/dt = dQ l jdt = 0, i.e. that the system be maximally 
integrable. the second one requires the system to be completely integrable, while 
the third one requires not only complete integrability, but also that there be stable 
equilibrium points. 

Thus, the HJ problem in this generalized form would solve not only a conjugacy 
problem, i.e. how to transform a given dynamical system into another one with 
a preassigned form, but one would find also the required transformation to be 
generated by a function S. 

/F rom a geometrical point of view, the first case requires the original dynamics 
to define a fundamental vector field of the natural foliation of the contact manifold 
after we have removed possible equilibrium points. In the second case the phase 
space is foliated by invariant cylinders, while in the third case the dynamics will 
preserve a foliation by tori. 

As already noticed, the conditions for the solution of the various problems are 
quite stringent, and one should expect that global solutions can be found only very 
rarely. Nevertheless, this more general perspective may be interesting because it 
could provide the Hamilton- Jacobi theory with a wider range of applications. For 
instance, it might be applied to the study of scattering problems, where now K 
would be the "comparison Hamiltonian" and H that of the system one is analyzing. 
It might be also applicable in Field Theory and General Relativity. In particular, 
one might consider applying it to Quantum Mechanics going beyond the usual 
WKBJ approximation. 

Remark 1. Let us observe that if we define, keeping the additional parameters 
fixed: 

dW 

Pi = -q-[ =Pi(q) (14) 
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and we form [14] the vector field: 

2(Q)3* = fJ(^(?))|l ("J 

it is immediate to show, using Eq.(5), that if: q = q (t) is an integral curve for 
the vector field X , then: (q (t) ,p (q (t))) is in turn an integral curve for Hamilton's 
canonical equations. Thus, Eq.(5) describes in a unified manner one particular 
family of the phase-space equations of motion. 

It is worth stressing that, already at this rather well-known level, the Hamilton- 
Jacobi theory establishes a deep connection between first-order partial differen- 
tial equations (PDE's) and systems of first-order ordinary differential equations 
(ODE's) [13]. 

1.2. The HJ Equation and the JWKB Method. Quite similar connections 
appear when one investigates the short- wavelength limit of "wave- like" equations 
(including also the Schrodinger equation) such as Hamiltonian optics as the short- 
wavelength regime of wave optics (the eikonal approximation [12, 54], but see also 
Ref. [36] for applications in Field Theory) and classical mechanics as the short- 
wavelength regime of wave mechanics [43] . 

In taking this limit, one starts with some differential operator of hyperbolic 
type on some manifold Q, passes through a Hamilton- Jacobi- type equation (the 
characteristics equation [35]) for a function S and, by substituting covectors (the 
"p"'s) for the first-order derivatives, arrives at a some Hamiltonian function on the 
cotangent bundle T*Q which yields also the dispersion relation of the wave motion. 
The associated Hamilton equations give rise, by projection on Q of the solutions, 
to the bi-characteristics [35] of the original differential system. Rephrased in an 
"optical" language, the characteristics describe the propagation of wave fronts, while 
the bi-characteristics describe that of rays. Note that "substituting covectors for 
first-order derivatives" , in the case in which the differential operator is homogeneous, 
is just a down-to-the -earth way of saying that we deal with what is known as the 
symbol [22] of the operator. 

Also, the JWKB method 1 [22, 23, 43, 45, 52] is widely used both in wave 
optics and wave mechanics to investigate several physical problems in the short- 
wavelength regime which, in wave mechanics, amounts to a leading-term expansion 
in powers of the Planck constant h. It is also closely related to the saddle-point 
approximation (plus one-loop corrections) in the path-integral approach to problems 
in Field Theory [2], Quantum [27] and Statistical Mechanics [31]. 

In this context, and concentrating for the sake of definiteness on the motion of 
a particle of mass to in a potential V (x) , one represents the wave function (in 
Gaussian form) as 2 : 

ip (x, t) = A (x, t) exp {iS (x, t) /h} (16) 

^his method was first introduced in the discussion of problems in wave propagation by lord 
Rayleigh back in 1912 and then applied to wave mechanics by H.Jeffreys in 1923 and, later on and 
simultaneously, by L.Brillouin, H.A.Kramers and G.Wentzel in 1926. 

2 Many authors [42] prefer to incorporate the prcfactor A into the definition of S, writing then: 
tp (x.f.) = exp(iS (x, t) /h, where the S of Eq.(16) is replaced by S — ihlnA and is no more real. 
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(with both A and S real). Substituting into the Schrodinger equation, one ends up 
with the (exact) coupled equations [22]: 

1 ,„„,2 t 7" / \ 55 ft 2 V 2 A 

— \VS \ +V(x) + — = — 17 

2m 1 1 v 1 dt 2m A y ' 

and: 

V(pv) + ^ = (18) 

where: p = p(x, t) = A 2 (x, t) and: v = v(x, t) = VS/m. Neglecting "quantum 
correction" proportional to h 2 on the r.h.s. of Eq.(17) yields the HJ equation (1). 
If S solves the time-dependent HJ equation, it is possible to show (see, e.g., App.4£? 
of Ref.[22] and Ref.[23]) that Eq.(18) is solved by 3 : 



P 



, / d 2 S 
det 



dx 2 dx 3 



(19) 



where the x^s are the initial coordinates, and hence the JWKB solution is: 



ip = y/pexp(iS/h) (20) 

As discussed in Ref.[22], Eq.(20) yields the JWKB approximation to the prop- 
agator (or Green function) for the Schrodinger operator, and the result becomes 
exact for quadratic Hamiltonians. Similar conclusions hold [2, 31, 33] in the path- 
integral formalism, where the saddle-point approximation 4 also becomes exact for 
quadratic Hamiltonians (or Lagrangians). 

1.3. A Geometrical Setting for the HJ Theory. A useful geometrical formu- 
lation of the HJ theory can be given as follows. Consider a symplectic manifold: 
M = T*Q with symplectic structure too = dpj Adq J . With any symplcctomorphism: 

0:M^M, 0*^0=^0 (21) 

(4>(<liP) ='■ (Qi-P)) we ma y associate a "graph" S^, i.e. a submanifold in the sym- 
plectic manifold M x M, the second factor being equipped with the symplectic 
structure — u>q: 

S = {{m,4> (m)) c M x M} (22) 

by requiring it to be Lagrangian w.r.t. the symplectic structure on M x M provided 
by: 

dpj A dq j - dPj A dQ 1 =: uj cj (23) 

By using the fact that: M x M = T*Q x T*Q ^T*(Qx Q), we can consider 
those Lagrangian submanifolds in T* (Q x Q) that can be written as graphs: 

dS :QxQ — >T*(QxQ) (24) 

When projects onto Q x Q we may consider it as the graph of a generating 
function S. 



3 The determinant on the r.h.s. of Eq.(19) is the well-known [27, 31] Pauli-Morette-Van Vleck 
determinant. 

4 In this case the Pauli-Morette-Van Vleck determinant is replaced by a functional determinant, 
whose definition requires careful regularization procedures [27, 56] that will not be discussed here. 
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In many cases M may admit of many alternative cotangent bundle structures, 
i.e. we may identify "alternative" submanifolds Q 1 such that: M = T*Q', with Q 
and Q' possible alternative "placements" of some external configuration space Q. 
A typical example is provided by: M = R 2 ™, where any identification of an affine 
subspace R n is a possible placement of an "abstract" M. n . 

^From this point of view, if we start with <fi and £^ we may look for a particular 
placement of Q in M that makes projectable (a submersion onto) Q x Q. 

Remark 2. The usual "colloquial" classification into "four possible sets of indepen- 
dent canonical variables" usually encountered in textbooks on Classical Mechanics, 
say (q, Q) , (q, P) , (Q,p) and (p, P) are exactly different identifications of the "con- 
figuration space" over which one constructs the cotangent bundle structure. There- 
fore we would have the corresponding symplectic structures: d (pjdqJ — PjdQ^ , 
d (pjdqi + Q^dPj) , d (PjdQ^ + q^dpj) and: d (q^dpj — Q^dPj) respectively. 

When the graph dS : QxQ — > T* (Q x Q) is viewed as a map: dS : QxQ T*Q, 
i.e. the second factor Q is considered as a family of "parameters" , we may obtain a 
regular foliation of T*Q. Then we may say that S is a "complete solution" of some 
associated H J equation, and in this case the Hamiltonian flow associated with the 
Hamiltonian K would preserve the foliation induced by 5" on T*Q. On the open 
dense submanifolds on which: dS : Q x Q — > T*Q provides a diffcomorphism we 
would have equivalence between the Cauchy problem in terms of initial data (qo,Po) 
and the boundary value problem in terms of (qo, Q). We notice that in this case: 

(dS)* c^o =: los = g^jW A d Q J (25) 

is a symplectic structure (cfr. Eq.(2)) on Q x Q, and hence it would allow for a 
Hamiltonian formulation on the space of " boundary data" rather than of the " initial 
conditions". The inverse image of the one-parameter group of evolution on T*Q 
would provide the "propagator" on the configuration space. 

By using the geometrical formulation of Quantum Mechanics a similar picture 
can be considered also for Quantum Mechanics in the Schrodinger picture. 

Let us restrict for convenience to a finite-dimensional Hilbert space IK. Again a 
transformation: 

$ : "K — ► 3i (26) 

will be associated with a graph: 

= {(tp, $ (VO) cMx 5C)} (27) 
If we consider onMxM the pseudo-Hcrmitian form: 

((^1,^2) I {<Pi,<Pi)) =■ (il>i\<Pi) ~ (ih\<Pi) (28) 

and noticing that, on £$, ipi^ = $ ("01,2), we find that S$ will be isotropic w.r.t. 
the pseudo-Hermitian form (28) iff $ is a unitary transformation. 

By considering the realification of JC, the Hcrmitian product decomposes into 
a real and an imaginary part, the former providing an Euclidean product and the 
latter a symplectic product. 
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As for the symplectic part, we have a situation similar to the one we considered 
previously, i.e. we may consider a generating function of the canonical transfor- 
mation and require afterwards that this transformation should preserve also the 
Euclidean product. In a somewhat simplified notation we may write then: 



Im (ip*dip — p*d<p) = dSq, 
and then we require, in addition: 

(dip\di/)) = (d(p\d(p) 



(29) 



(30) 



The resulting transformation would be a Kahlerian transformation derived from 
the generating function For: ip = Utjj, tp* = ip*U*, S((p*,ip) = —<p*Uip: 



d(p* A d(p- di'* A dip 
For example, a Hadamard gate: 



-d (tpdtp* + ip*dip) 



U = U » = T2 



1 1 
-1 1 



has the generating function: 



(31) 



(32) 



(33) 



1.4. Quantum HJ Equations. In the Heisenberg picture, one may devise a di- 
rect approach to the HJ problem by replacing classical variables with Hermitian 
operators. Indeed, shortly after the advent of Quantum Mechanics, various efforts 
were made [19, 21] to formulate both the theory of Canonical Transformations and 
hence also the Action Principle [50] and the Hamilton- Jacobi equation in operator 
terms from the very beginning . 

Following, e.g., the scheme of Eqs.(l) to (3), one might be tempted to "promote" 
these equations to operator equations defining a (quantum) canonical transforma- 
tion as: 



P, 



d_ 

dq 1 



(q,Q, 



P, 



_d_ 



S 



(q, Q, t) 



and a (quantum) HJ equation as: 

^ dS 
H\q,—,t 



t ) = 



(34) 



(35) 



Due to operator-ordering problems, these equations are obviously ambiguous. 
According to Jordan [29, 30] and Dirac [19, 20] the ambiguity should be resolved 
by requiring S and the operators in Eqs.(34) and (35) to be "well ordered", i.e. 
by requiring all the "uppercase" operators (the Q's) to stay to the right of the 
"lowercase" ones (the q's). This implies that the "generating operator" S should 
be of the general form [46] : 



S 



(q, Q, t) 



= ^2 fa (q, t) 9c 



(Q,t) 



(36) 



5 See Ref. [50] for a systematic exposition of this approach. 
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for suitable functions f a and g a . The "well-ordering" procedure must be applied 
also, using when necessary the commutation relations, to the Hamiltonian operator 
H in Eq.(35). 

The authors in Ref.[46] have devised a procedure for converting the operator 
equation (35) into a c-number equation that we will illustrate here on a simple 
example 6 , namely that of a (non-relativistic) ID particle of mass m subject to a 
scalar potential V (q). The Hamiltonian is then: 

^2 



and the (operator) H J equation becomes (S = S(q,Q,t)): 

2 



1 

2m 



OS 



(37) 



(38) 



Sandwiching between eigenstates \q) and \Q) of q and Q respectively, the "well 
ordering" prescription leads to: 

(39) 



(q\S(q,Q,t) \Q)=S{q,Q,t) (q\Q) 

i.e. to a (uniquely defined and not necessarily real) c-number function S(q,Q,t). 

Sandwiching between the same eigenstates the quadratic term on the l.h.s. of 
Eq.(38) is a bit more complicated. Explicitly: 



OS 



(40) 



Now, the standard canonical commutation relations imply, for any function G = 
G{q) : 

[G,p\=iH— (41) 
Using then the first of Eqs.(34) (i.e.: dS/dq = p) one finds [46]: 
dS dfp = dfp dS ih d 2 fp = dfp (q, t) df a (q, t) ~ . n &ffi(M,t) (42) 



dq 2 



which is again a "well ordered" expression, and hence: 



dS 



\Q) 



S(q,Q,t) 



— 

dq- 



S(q,Q,t) 



(q\Q) 



(43) 



Dropping then the common factor (q\Q) one obtains the c-number equation 8 : 



1 

2m 



d_ 

dq 



-ih^S(q,Q,t) 



V(q) + ^- t S(q,Q,t)=0 (44) 



which is completely equivalent to the operator equation (38). 



6 Referring to Ref. [46] for a more general discussion. 
7 J dq \q){q\ = I, {q\q') = S (q — q'), and similarly for Q. 

8 As Eq.(44) (as well as the classical equation (1)) contains only derivatives of S, any solution 
will be ambiguous by the addition of a constant term. This ambiguity, which is totally irrelevant 
at the classical level, will turn out instead to be useful in what follows. 
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Eq.(44) is precisely the equation that results [26, 46] from the Schrodingcr equa- 
tion (in the variables (<?,£)) by expressing the wave function ip as (see footnote 
2): 

il){q,Q,t)=expUs{q,Q,t)\ (45) 

It provides then a solution of the Schrodingcr equation depending (in an essential 
way) from the additional parameter Q, just as the propagator K(q,Q,t) [22, 49], 
which obeys the same equation, docs. The propagator is known [22, 31, 49] to obey 
the boundary condition: K (q, Q,t) — > S (q — Q) as t — > and, in order to complete 
the identification, one has to check (the equation being first-order in time) that 
there is a solution of the form (45) which does the same. For example, from the 
well-known 9 result [31, 49] for the propagator of the ID harmonic oscillator with 
mass m, proper frequency u and Hamiltonian: 

H=^- + \mu 2 q 2 (46) 
2m 2 

one can check directly that: 

mu T , , , . „ ^ ih, ( 2wih sin (ut)\ .,„. 

S = ^ + Q ) C ° S ™ - 2 ^ + ~2 ^ ( TflLO ) ^ 

does indeed solve Eq.(44) with the appropriate boundary condition 10 . It is then 
immediate to see that the "well-ordered" solution of the operator equation (36) will 
be: 

S = - . . [f + Q 2 )cos(cot)-2qQ + — In I 48 

2sm(wt) LV / J 2 \ muj J 

Remark 3. When substituting from Eq.(48) into Eq.(36) we need to square the 
derivative: dS/dq = muj(qcos(ujt) — Q)/sin(uit). This brings about terms that 
are not " well-ordered" . Bringing them in the correct order [46] requires using the 
(exact) commutation relation: q , Q = —ihsm(ujt) /mu. The time derivative of 
the last term in Eq.(48) (see also footnote 8) does then the job of compensating for 
this operation. 

The short-time limit of Eq.(48) is: 

S t ~0 2t( q +Q - 2qQ ) + ^{—) (49) 

which is the same as for the free particle. The harmonic potential does not con- 
tribute to this short-time limit. As remarked in Ref. [46], the same will happen more 
generally for any nonsingular potential V(q). So, the result (49) will have a more 
general significance, going slightly beyond the case of quadratic potentials. 



9 As for any quadratic Hamiltonian. 

10 In the limit uj — > 0, insertion of Eq.(47) into Eq.(45) reproduces of course the kernel for the 
ID free particle. 
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1.5. Comments and Plan of the Paper. Going back now to our general dis- 
cussion, the attempt to recover the original PDE from the equations of motion of 
the "rays", i.e. Hamilton's equations, is what is usually called the "quantization 
problem" . 

This interplay between PDE's, associated first-order PDE's (HJ-type equa- 
tions) describing the propagation of " wave- fronts" and the corresponding, "ray" 
ODE's (Hamilton's equations) has been widely investigated in various branches of 
theoretical Physics. Moreover, because of the Hamiltonian-Lagrangian correspon- 
dence, also the calculus of variations appears in this interplay. 

^From all these remarks it should be clear that the H J theory is very rich in an- 
alytic and geometric ideas, and that it unifies apparently diverse topics like higher- 
order PDE's, first-order PDE's, ODE's and the calculus of variations. 

In this paper we will try and use some of our experience with Quantum Mechanics 
to formulate and give a geometric presentation of many problems which, born in 
a quantum setting, are of more general validity in the framework of the Hamilton- 
Jacobi theory. Besides many original papers on this vast subject, we shall rely on 
some work by A.Vinogradov [53], some more recent work by Grabowski and Poncin 
[25], a previous paper of ours [38] and recent paper on the Hamilton- Jacobi theory 
in a Lagrangian setting [14]. We shall discuss the following topics: 

-Reviewing briefly how the H J problem can be formulated in geometric terms 
on the cotangent bundle T*Q of a smooth manifold Q, we shall discuss how the 
same problem can be formulated on the tangent bundle TQ, hence in a Lagrangian 
setting. 

-How one can pose a generalized H J problem for differential operators of any 
order, giving another coordinate-free characterization of the HJ equation. 

-We shall try and discuss to which extent, instead of differential operators acting 
on functions, one can consider differential operators acting on sections of vector 
bundles, thereby obtaining "wave-like" equations that are not scalar, like the Pauli 
and the Dirac equations. 

-Just as in Quantum Mechanics one poses a joint eigenvalue problem for two (or 
more) observables, we shall discuss how one can pose a "joint HJ problem" for 
more than two functions on the cotangent bundle, and, finally, 

-After "revisiting" briefly the geometrical formulation of the timc-dependent 
Hamilton- Jacobi theory, whose proper setting [38] is on the cotangent bundle T* (Q x 
we shall discuss how one can obtain generalizations thereof when the action of the 
Abclian group R is replaced with that of a general Lie group G (a "HJ problem on 
a Lie group", then). 

2. A Geometrical Setting for the Time-Independent Hamilton-Jacobi 
Theory on the Cotangent and on the Tangent Bundles. 

2.1. Preliminaries. Let us begin by recalling [38] some preliminary notions. The 
cotangent bundle T*Q carries with it the canonical (exact) symplectic structure: 

uj a = d6 (50) 
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where, in local coordinates: 9q = pidq 1 , and hence: ujq = dpi A dq l . 

Let then a G X* (Q) be a one-form on the base manifold Q. Again, in coordinates: 

a{q) = a i (q)dq i (51) 

With the one-form a we can associate the map (the graph of a): 

<p a :Q^T*Q by: Q 3 q » {q,a (q)) g T*Q (52) 

The image of a, T [a] will be defined as: T [a] = ip a (Q). Of course: dimT [a] = 
dimQ = n. It is clear from Eq.(52) that: 

ir o ip a = Id Q (53) 

where: tt : T*Q — > Q is the canonical projection. This shows that ip a gives a global 
section of T*Q. T [a] will be therefore an n-dimcnsional transversal [38] submanifold 
of T*Q and ip a will be an embedding [41] of Q into T*Q. Also, a can be recovered 
from the canonical one-form 9q via the pull-back: 

tp* a 6 = a (54) 

We can consider also the map: 

V Q =: <p a o 7T : T*Q -> T*Q fey : T*Q 3 (g,p) h-> (q,a(q)) G T*Q (55) 

which is a base- invariant translation along the fibers that, for every q £ Q, "shrinks" 
the whole fibre T*Q to the point a (q). Then, using Eq.(54), one obtains at once: 

■0*^0 = n*da (56) 

and we conclude [38] that the graph of a, besides being always transversal to the 
fibers, will be also a Lagrangian submanifold [38, 41] of T*Q if and only if a 
is closed. As discussed in Rcf. [38], the converse is not true, i.e. a transversal 
Lagrangian submanifold of T*Q which projects down to Q under the canonical 
projection need not be the graph of a closed one-form. 

Being closed, a will be locally exact, i.e.: a = dW for some function W £ 3 (Q), 
at least locally. Whether or not such a function exists globally will depend on 
whether or not H 1 (Q), the first de Rham cohomology group of Q, is trivial. 

2.2. The HJ Theory on the Cotangent Bundle. After these preliminaries, 
let us re-consider now the time- independent HJ equation, Eq. (5), for Bamilton's 
characteristic function, i.e.: 

l(r,T)-B (57) 



dq 

W can be either a particular solution of the HJ equation or a complete in- 
tegral (with all the possible intermediate cases in between). In the latter case: 
W = W (q; a) depending (as already discussed, in an essential way) on n additional 
variables collectively denoted as a. Whenever necessary, we will denote as W a (q) 
the function that obtains by keeping the a's constant. Bence: W a £ 5 (Q) Va. 

It is clear that, considering the graph of the exact one- form dW a , the image of 
dW a is a Lagrangian submanifold and Eq.(57) can be rewritten as [38]: 

dW* (H-E) = (58) 
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or, equivalently, as: 

dW* (dH) = (59) 

As it stands, Eq.(58) looks just like a different way of rewriting the standard 
H J equation in a different language. However it allows for a deeper geometrical 
interpretation of the HJ theory, and allows also for some interesting generalizations. 

First of all, the image r a of dW a is a regular submanifold [41] of T*Q. As such, 
it can be described, locally at least, as the zero-level set of n independent functions 
ha,-,fna € ${T*Q), i.e. such that: 

dfla A ... Adfna ^ 

T a = {m=(q,p)eT*Q\ f ja (m),j = 0, l,...,n} 1 ' 

and: 

{fia,f ja }(m) = (61) 

where {., .} is the Poisson bracket in the space of functions on T*Q associated with 
the symplcctic form loq. 

Eq.(61) implies that the Hamiltonian vector fields Xj associated with the f/s, i.e.: 

ix ja uo = -df ja , j = 1, ...,n (62) 

are all tangent to the submanifold (indeed [38] they span the tangent space T m T a 
at each m € T a ). Moreover, denoting by Xh the Hamiltonian vector field that 
describes the dynamics, we obtain, contracting both sides of Eq.(59) with the Xj a 's: 

= [ix ia (dW: (dH))] (m) = [i Xja dH] (m) = 
= (L Xja H) (m) = - (L XH f ja ) (m) = {H, f ja } (m) 1 ' 

which proves that Xh is tangent to r a Va. Of course Xh is also tangent to the 
(2n — l)-dimcnsional energy surface 11 : 

S B = (H - E)" 1 (0) (64) 

If dH ^ 12 , Yje is a regular submanifold that we can assume without loss of 
generality to be also connected. If it is not, we can always restrict the discussion to 
each connected component separately. Each T a is contained in the energy surface 
(64) and, in general [38], they will provide a (n — l)-parameter foliation of the 
energy surface, with the dynamical vector field Xh being tangent to all the leaves 
of the foliation. This is basically the geometrical content of the fact that, besides 
the energy E, a complete integral depends in an essential way on n — 1 additional 
parameters. 

This formulation of the HJ problem suffers however of some limitations, as the 
following example shows. 

Example 1. Consider: Q = S 1 with coordinate q, < q < 2ir. The cotangent 
bundle can be given coordinates (g, p) , with pgK, and can be viewed as a cylinder. 
In this case: 9o = pdq and: uo = dp A dq are both well-defined. The vector field: 
X = pd/dq is Hamiltonian with the Hamiltonian: H = p 2 /2, and the "energy 



11 Remember that, in the time-independent case, the energy E must be included among the pa- 
rameters on which a complete integral depends, leaving actually only n— 1 independent parameters. 

12 Critical points of H, being invariant sets for the dynamics, can be handled separately. 
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surfaces" are pairs of circles on the cylinder: = yq, ±\/2Ej , E > 0. The 

associated HJ equation: (dS/dq) 2 = E has of course no global solutions for E > 0, 
but the energy surfaces are nonetheless the graphs of the closed but not exact one- 
forms: otE,± = ±v2Edq, and the equation: (p* a (H — E) = for an unknown closed 
one-form a is globally defined and has precisely cxe,± as solutions. 

This suggests that one relaxes the requirement that the pull-back of H — E in 
Eq.(58) be via an exact one- form, replacing it with the weaker request that it be 
via a closed but not necessarily exact one-form, i.e. that one replaces Eq.(58) with: 

a* (H — E) = 0, a e X* (Q) , da = (65) 

The graph of any such form will be of course again a transversal Lagrangian sub- 
manifold, and we can re- formulate the "geometric HJ problem" as follows: 

• A solution of the HJ equation for a given Hamiltonian and a given energy is 
a transversal Lagrangian submanifold within the energy surface, obtained as 
the graph of a closed one-form on Q, and 

• A complete integral is a foliation of T*Q by such solutions for all physically 
accessible values of the energy, which implies also a foliation of each energy 
surface as well. 

Note that (cfr. Eq.(54)), as: a*cu = a*dd = d(a*8 ) = da, Eq.(65) can be 
replaced by the equivalent one: 

a* (H - E) = 0, a e X* (Q) , a*uj = (66) 

If, as a further step, one gives up the requirement of transversality, one can pose a 
V HJ problem" (no more a "HJ equation") consisting in the search for foliations of 
each energy surface simply by Lagrangian submanifolds (that, being not necessarily 
transversal, can exhibit caustics [4]). This has the advantage that the full set of 
canonical symmetries can be implemented as symmetries of the HJ problem (i.e. 
maps that map solutions into solutions). We will not insist on this point, but refer 
rather to the literature [38] for a more complete discussion. 

2.3. The HJ Theory on the Tangent Bundle. We turn now briefly to the 
Lagrangian context. The relevant carrier space is now the tangent bundle TQ, 
with local coordinates (<f , m 4 ) , i = 1, n, which carries no pre-assigned symplectic 
structure but, if a Lagrangian £ € J (TQ) (assumed here to be regular 13 ) is given, 
can be endowed with a Lagrangian symplectic structure loh defined by: [1, 3, 41]: 

BL 

uc=d9i; 6 L = —dq l (67) 
ou 1 

again, for simplicity, in local coordinates. The Euler-Lagrange equations can be put 
in "Hamiltonian" form as: 

ir L uja = -dE L (68) 

where L^ is the second-order [41] vector field describing the dynamics, the "energy 
function" E^ is given by: 

E C = (L A -1)L (69) 



We refer to Ref. [41] for a discussion of the case of singular Lagrangians. 
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and, finally, A is the dilation (Liouville) field along the fibers: 

A = ^ ™ 
Sections (actually, global sections) of the tangent bundle are now provided by 
vector fields on the base manifold, just as one-forms did the same job for the cotan- 
gent bundle. If X € X (Q) is any such vector field, given in local coordinates by: 
X = X l (q) d/dq\ then X will define the map 14 : 

X:Q^TQby: ^ (q\ X* (q)) (71) 

satisfying: 

7T o X = Id Q (72) 

where: n : TQ — > Q is the canonical projection, and hence X is a section of TQ. 

With reference to Eq.(66) we can then define a Lagrangian HJ problem as the 
search for all vector fields X £ X (Q) such that: 

X* {E L -E)=0 and : X*uj l = (73) 

As such, the HJ problem on the tangent bundle can be viewed as the search, 
instead of a single function as in the case of T*Q, of a " vector- valued function" on 
Q. 

As a (Lagrangian) counterpart of the Remark that was made in Sect.l, we have 
now the following [14] 

Remark 4. If: X g X (Q) , X : Q -> TQ is a solution of Eq.(73) and 7 : R -> Q is 
an integral curve of X, i.e.: 7 = X 07, then: 7 : K — > TQ solves the Eulcr-Lagrangc 
equations for the Lagrangian £, i.e.: 



107 = 7 ==> L,c 07 = ^07 (74) 

The converse statement, however (i.e. if the integral curves 7 of a vector field 
X g X (Q) are such that 7 are integral curves for the Lagrangian vector field 
for a given Lagrangian L, then X is a solution of Eq.(73)) need not be true, as the 
following example [14] shows. 

Example 2. The dynamics of the free particle in 
regular Lagrangian: 



1 r 



L =2 



K) 2 + KT 



2 can be described by the 
(75) 



with associated geometrical objects: 



Tr = U 



1± 

dq 1 



and: 



= ^dq 1 + u 2 dq 2 ; Ll>c = du 1 A dq 1 + du 2 A dq 2 



(76) 



(77) 



14 With some abuse of notation, we are using here too the same symbol to denote the vector 
field and the associated map. 
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The two-parameter family of vector fields: 

, d kq 2 - 1 d 

X = k w + —^ (78) 

satisfies the assumptions of the above (putative) converse statement, but: 

kn 2 — I 

X*(^) = -^—fdq 1 Adq 2 ^0 (79) 

*d(X*(E z ))?0 (80) 



and: 

'kq 2 - I 



'I 1 



Hence, both conditions of Eq.(73) will be violated. 

Guided by this example we will stick to the definition of the Lagrangian Hamilton- 
Jacobi problem as defined by Eqs.(73), i.e. as the search for all vector fields X £ 
X (Q) such that (using again the same notation for the associated maps Q — > TQ): 

X*uj L =d(X*E L ) = (81) 

The first of these equations implies, of course, that ImX be a Lagrangian sub- 
manifold of TQ. Furthermore, as: = X*loh = X*(d6 Cj ) = d(X*dc), every point 
has an open neighborhood U C Q where there is a function W £ £F (U) such that: 
X*6 C = dW in U. 

We have already shown in a previous Remark that if X is a solution of the 
Lagrangian HJ problem, then it satisfies Eq.(74) (while the converse is not true) 15 . 
It has been proved in Ref.[14] that X is a solution of Eq.(73) iff the following 
diagram: 

TQ ^ T (TQ) 

X t T TX (82) 

Q — ► TQ 
x 

commutes, i.e. iff: 

T L oX = TXoX (83) 

Eq.(83) is a PDE for the unknown vector field, or "vector- valued function", X 
which replaces the PDE for the scalar function W. Once a solution is found, one 
has however to check whether or not it satisfies the conditions (81) as well, in order 
for it to be a genuine solution of the Lagrangian HJ problem. 

It may be useful to derive the expression for the PDE (83) in local coordinates. 
We can write X and Tc as: 

X = X 1 -£t; X 1 =X i (q) , 

" 9 (84) 



r L = 

where: 



dL d 2 L 



15 In Ref. [14] the problem of finding the solutions of Eq. (74) is also called the generalized 
Lagrangian HJ problem. 
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and H %3 is the inverse of the Hessian matrix: 



d 2 L 

A direct calculation shows that 



Hij - du^J (86) 



(ftvi \ ft 

w x3 ~ at{q ' x) )d^ (87) 

This is a vertical vector field along X whose vanishing implies: 



-JP -a t (q,X) = O i i=l,...,n 
and this is the required local form of the PDE (83). 

One of the main results of Ref. [14] can now be rephrased as follows: 
Proposition. The following statements are equivalent: 

1. X is a solution of the Lagrangian HJ problem. 

2. Besides being a Lagrangian submanifold of TQ, hnX is also invariant under 
the dynamics represented by i.e. T& is everywhere tangent to ImX. 

3. The integral curves of with initial conditions on ImX project onto the 
integral curves of X. 

In the context of the tangent bundle, the notion of a complete solution of the 
Lagrangian HJ problem can be posed as follows: 

Definition. A complete solution of the Lagrangian Hamilton- J acobi problem is 
provided by a family of solutions {X\} X£A , A an open set in R n , such that the map: 

$:QxA^TQ by: $ (q, A) = X\ (q) (89) 

is a local diffeomorphism. 

It follows from the definition that a complete solution yields a foliation 16 of TQ 
with leaves transversal to the fibers, and that the Lagrangian vector field is 
tangent to the fibers of the foliation. 

Example 3. Consider the two-dimensional harmonic oscillator with the stan- 
dard Lagrangian: 



£ = i 
2 



The dynamical vector field is: 



(90) 



L dq 1 dq 2 ^ du 1 ^ du 2 ^ 



and: 



uj l = du 1 A dq 1 + du 2 A dq 2 



(92) 



'This is ensured by the request that <E> (q, A) be a (local) diffeomorphism. 
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It is known that the functions: 



fo = u l u 2 + q x q 2 

A = (- 1 ) 2 + (? 1 ) 2 m) 
/ 2 = (" 2 ) 2 + (^ 2 ) 2 ( } 

f 3 = q x u 2 - q 2 u x 

are all constants of the motion, not functionally independent, of course. Let, e.g., 
fo = C, f\ = 2E\, /a = 2E 2 , /3 = I. In particular, f\ and f 2 are in involution: 
{/i;/2} = 0. The PDE equations (88) read now: 

w X' = -J,i=l,2 (94) 
and it is easy to check that the four two-parameter families of vector fields 17 : 

X El , E2 = ±^2 El (gi) 2 A ± ^2E 2 -(q 2 ) 2 ^ (95) 
do indeed solve them. As: 



(X ElE2 )* uc = ±d^/2E 1 - {q 1 )' A dq 1 ± d\J2E 2 - (q 2 f A dq 2 = (96) 

and: 

(Xe 1 ,e 2 )* dEc = 



>2E X - {q 1 Yd^2E 1 - (qiy + ^2E 2 - {q 2 Yd^2E 2 - {q 2 Y + q l dq l + q 2 dq 2 = 

(97) 

each family (95) is indeed a complete solution of the Lagrangian HJ problem for 
the two-dimensional harmonic oscillator with the standard Lagrangian (90) . 

Remark 5. It is useful to stress here that the PDE (88) (or, for that matter, 
Eq.(94) in the case of the harmonic oscillator), which defines [14] the "generalized" 
Lagrangian HJ problem, depends on the Lagrangian only through the vector field 
Tjr,, while the HJ problem that we are discussing here depends also on additional 
structures derived from the Lagrangian, namely the symplectic form and the 
" energy function" Ec . 

It is known [44] that the so-called "Inverse Problem in the Calculus of Variations" , 
i.e. the problem of whether or not a dynamics described by a given second-order 
vector field on TQ admits of a Lagrangian description can have no solutions at 
all, only one solution 18 , or more than one solution, leading then to (genuinely) 
alternative Lagrangian descriptions for a given dynamics. 

If this is the case, vector fields that are solutions of the Lagrangian HJ problem 
(or families thereof, yielding complete solutions) for a given Lagrangian need not 
be such for an alternative Lagrangian description, while remaining, according to 
what has just been said, solutions (or complete solutions) of the generalized HJ 
problem, as the following example shows. 



17 It is easy to recognize that the components of the vector fields are obtained by expressing 
2 in terms of the coordinates on the base manifold and of the two parameters -Ei : 2- 
18 Apart from the addition [44] of essentially trivial "gauge terms". 
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Example 4. Consider again the two-dimensional harmonic oscillator, but now 
with the alternative Lagrangian [44]: 



£.1 



1 



K) 2 -( u 2 ) 2 -(/) 2 + ( g 2 ) : 



The associated structures will be now: 



ujh 1 = du 1 A dq 1 — du 2 A dq 2 



(98) 



(99) 



while, of course, the dynamical vector field and the PDE (94) will be the same. 

The vector fields (95) will satisfy again (Xe 1 e 2 )* = as well as: (Xe 1 ,e 2 )* dEn 1 
and hence the Xe 1 e 2 's will provide a complete solution of the HJ problem for the 
Lagrangian (98) as well. 

If one considers instead the Lagrangian [44] : 



we have: 

and one finds, e.g.: 
(Xe 1 e 2 )* ^2 = 



± 



L-2 = u l u 2 — q 1 q 2 



E L2 = u 1 u 2 + q 1 q 2 
ujh 2 = du 2 A dq 1 + du 1 A dq 2 



± 



2E 2 - {q 2 



2E 1 - ( q iy 



(100) 
(101) 

dq 1 A dq V (102) 



(as well as: (Xe 1 e 2 ) (dEj^ 2 ) ^ 0). Hence, the Xe 1 e 2 's will be no more solutions 
of the HJ problem for the alternative Lagrangian (100). 



3. The Generalized Hamilton-Jacobi Problem for Differential Operators. 

3.1. Differential Operators and Principal Symbols. In order to deal with 
differential operators on manifolds, we will first review differential operators on R". 
We will do this by providing an algebraic characterization that will allow us to deal 
with differential operators on arbitrary manifolds. 

We consider then the algebra A — 3 (W 1 ) of smooth functions on R™. A differ- 
ential operator of degree at most k is defined as a linear map: L>( fe ) : A -> A of the 
form: 

D [k) = E-^f^; 9a eA (103) 

|cr|<fe a 

where we have introduced multi- indices: a = (ji, i n ) , |er| = i\ + i 2 + ... + i n , 
and: 

— = : : (104) 

dx a dx\ 1 dx 2 2 ...dx z n" v ; 

It is possible to give an algebraic characterization, appropriate to arbitrary mani- 
folds, in the following way. With functions / E A we associate differential operators 
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/ of order zero that act by multiplication, i.e.: fg =: fg on all smooth functions. 
We notice that the commutator bracket gives: 



_d_ 

dxi 



and, more generally: 



f 



E 



df_ 

dxi 



dx r dx v 

for |r| > (strictly) and some set of constants c T . It follows then easily that 

'Si 



(105) 
(106) 



dx a 



(107) 



\a\<k 

is a differential operator of degree at most k — 1. Iterating the procedure for a set 
of k + 1 functions fo, fi, fk, we find: 



) — i fk 







(108) 



The converse statement holds also true, namely, a linear operator which does not 
increase the support and satisfying the property (108) on any set of k+1 elements in 
A will be a differential operator [5] . By means of this algebraic characterization it is 
then possible to define differential operators as linear maps satisfying the property 
(108). 

We notice that: 



D {k \D^ = D {k+j - 1] , k+j > 1 



(109) 



• Setting: k = j = 1 we obtain that differential operators of degree at most one 
are a subalgebra and that A, as an Abelian subalgebra of operators of degree 
zero, is an invariant subalgebra thereof. 

• Differential operators of degree at most one are derivations of A if they are 
zero on constants. 

• Derivations arc a subalgebra (the subalgebra of homogeneous differential op- 
erators of degree one) . 

Remark 6. Considering the algebra A of smooth functions and the algebra DerA of 
the derivations on A, and noticing that: [X\ + fi,X2 + $2} = [^1,^2] + ixi/2 — 
Lx 2 fi, X\,X2 € DerA, /i,/2 G A, we can form a semi-direct product of Lie 
algebras by setting: 



Pa,/i) , (X 2 , fi)] = ([X U X 2 ] , L Xl f 2 - Lxji) 



(110) 



The enveloping algebra of this Lie algebra will be isomorphic with the algebra of 
differential operators. 



If D is a differential operator, then for any two functions / and g: 



fdd,D] 



f,D 



(111) 



It follows then that: 
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Proposition. If D^ k > is a differential operator of degree at most k : the expres- 
sion: 



D« /i 



■fk 



(112) 



is a function which is symmetric with respect to all the permutations of fi, ...,/&.■ 
We can set up an equivalence relation " ~ " among differential operators of the 



same degree, say fc, by saying that D^ ~ D^' iff: 



,(fe) 







■■■fk 








■■■fk 



V/i, ■ •■,/& 



(113) 



(114) 

(115) 
(£>(«)) = 

(116) 



The equivalence class is what is called the principal symbol of the differential op- 
erator. The set of principal symbols is a commutative algebra, this following from 
the fact that two operators of degree k are in the same equivalence class iff their 
difference is a differential operator of degree < k. The set of the symbols of the 
differential operators of degree k will be denoted as (Q). The action on a set 
f = (/i, .., fk) of function of the principal symbol of an operator of order k 
will be denoted as a p (D^)) (f) and, of course: 

a p (Z#>)(f)=[...[£« /i] ...fk 

We note that: 

(a p (d«) a p (f) = o p (d« • DM) (f) 

with, now: f = (/i, fk+ m ), which shows also that: a p (D^ k n a p 
a P (£>('»)) a P (£)(*)), as well as that: 

[Di,Da-D 3 ] = [D 1 ,D 2 ]-D 3 + D 2 -[D 1 ,D 3 ] 

for any three differential operators D\, D 2 and D 3 . 

It is possible to define a Lie algebra product on principal symbols by associating 
with any two of them the principal symbol of their commutator. If we denote by: 
a p : Dif '/W (Q) — > (Q) the map that associates a symbol with an operator, 
we can define: 

{a p (Z?0 , a p (U a )} =: <r p ([£>!, D 2 }) (117) 

In this way we define a Poisson bracket on the commutative algebra of the principal 
symbols of differential operators. Moreover, we have: 

{a p (DO , a p (D a ) a p (D 3 )} = a p ([Di,D 2 ]) a p (D 3 ) + a p (D 2 ) a p ([D 1} D 3 ]) 

(118) 

which shows that {<r p (D) , •} is a derivation on the commutative algebra of the 
principal symbols. 

It is not difficult to see from Eq.(114) that one can define (and in an unique way) 
a symmetric contravariant tensor of rank k via: 

o p (DM) (f) = (dh,df 2 ,...,dfk) (119) 



(120) 



;dfk) 

evaluated on the symmetrized product df\ ® df 2 ® ... ® dfk- 

For example, for a homogeneous second-order operator of the form: 



2 , 



22 



G. MARMO G. MORANDI AND N. MUKUNDA 



z 3% ) we find: 

2 



This characterization of principal symbols by means of symmetric contravariant 
tensor fields allows us to conclude that the Poisson bracket on principal symbols is 
isomorphic with the Lie algebra product on symmetric contravariant tensors defined 
by the Schouten bracket [48]. 

Remark 7. If Qis parallclizable, there will be a global basis, say {X±, ...,X n ) of 
vector fields, and we can consider the 7 (<5)-module of contravariant tensor fields 
generated by them. For example, the monomial X^ <g> ... <g> Xi k defines a differential 
operator of order k by setting: 

D^(f) = L Xii (L Xta ...(L Xi j)) (122) 



which corresponds to the differential operator 



, which is of order 



k commutators 

zero. 



A further identification is possible by considering the principal symbol a p (Z)( fe )) 
as a fiberwise polynomial function f D (k) on T*Q defined as: 

f DW =:DW(flo®flo®...®flo) (123) 



For example, with: 9q =Pidq l , the rank-two tensor (121) leads to: 

f D m = l alJ P l Pj ( 124 ) 



Let us stress the fact that what Eqs.(119) and (123) show is that it is possible to 
characterize principal symbols completely in tensorial terms. 

Remark 8. We are committing a slight abuse of notation in making this definition, 
as wc arc contracting tensor fields on Q with forms on T*Q. However, the abuse is 
justified by the fact that 9q is a semi-basic one-form. 

With this association, the Poisson bracket (associated with d6o) on polynomial 
functions defined by contracting D^ k > with covariant tensor fields defined by means 
of powers of 9 turns out to be isomorphic with the "abstract" Poisson bracket (117) 
defined on principal symbols of differential operators. 

We should stress here that this realization in terms of functions that are poly- 
nomials in the momenta depends on the semi-basic one-form we have used (even 
though 9q is a natural one- form [41] on T*Q). It would be possible to consider other 
semi-basic one-forms whose exterior derivative would be a symplectic structure. 
For example, we might consider, for any non-singular numerical matrix K — \\K l A\: 



9 K = Pi K l jdq j 



(125) 
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whereby: d9x = dpiK 1 j A dq 3 will be a symplectic structure. Then, again using 
the rank-two tensor (121), one finds the quadratic form: 



D (2) {9 K 9) ' 



} k) = -Pi 



Ka K) 



Pj 



(126) 



1 = K % j^j. In this way, again the associated Poisson algebra would be 
isomorphic with the abstract Poisson algebra defined by means of the principal 
symbols. These alternatives may turn out to be relevant if we would like to consider 
bi-Hamiltonian systems. In such a situation, semi-basic one-forms 9 such that: 

L v 9 = dF (127) 

which implies: 

i r dO = -dH; H =: i T 9 - F (128) 
would be of particular interest. 

Example 5. On T*M. 2 with coordinates (q 1 ,q 2 ,Pi,P2) the dynamics of the 2D 
isotropic harmonic oscillator can be represented by the vector field: 



_ d d , <9 2 d 

r = P!t~T + P' 2 ^2 ~ 1 a 1 a - 

aq L oq^ opi op2 



(129) 



which is Hamiltonian w.r.t. the "canonical" symplectic form uq = d9o = dpi A 
dq 1 with the "standard" Hamiltonian: H = ^.^{pi) 2 + (g l ) 2 )/2. It is also Hamil- 
tonian (and hence bi-Hamiltonian) w.r.t. (among others [44] ^ the symplectic form: 

d9; 9=p 2 dq 1 + Pidq 2 (130) 



which is of the form (125) with: 



K 



id where: 



L T 9 = dF; F = p lP2 



(131) 
(132) 



3.2. Hamilton-Jacobi-Type Equations Associated with Symbols. Up to 

now we have seen that with the principal symbol of any differential operator we 
can associate both a symmetric contravariant tensor field and a polynomial func- 
tion on T*Q. Clearly, if we consider this polynomial function fo, it is possible to 
construct a first-order PDF by setting, as in Sect. 2: 

(dS)*(f D ) = (133) 

or, more generally: 

(dS)* ( f D ) = c (134) 

thus defining a PDE of the Hamilton- Jacobi type. 
^From what we have said, we may set, equivalcntly: 



5' 



(135) 



so that we can define the HJ equation associated with D^ k > directly on the config- 
uration manifold Q, without making recourse to the cotangent bundle. 
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Out of the cotangent bundle representation of the principal symbol a p (D), i.e. 
by means of the polynomial functions fn on T*Q, it is possible to associate also 
Hamilton's equations constructed out of fo, namely: 

i r d9 = -df D = -d (D (0 O ®9 ® ... ® 6> )) (136) 

In this way, by means of the principal symbols, we are able to construct a first- 
order PDE of the Hamilton-Jacobi type as well as Hamilton's equations that define 
the vector field T. Let us look now at a couple of examples. 

Example 6. The Schrddinger operator. The operator associated with the evolu- 
tion equation is: 

Ds = ih^- + ^X7 2 -V(r) (137) 
ot 2m 

The associated HJ equation may be written in the form: 



D S ,S 



r, 2 



.s 



E (138) 



2m VSf=E (139) 

We find that the principal symbol, or the HJ equation associated with it, contains 
no information on the potential nor on the time evolution. 

Example 7. The Klein- Gordon operator. A similar situation prevails for the 
Klein- Gordon operator: 



i an. 



D KG = V 2 --^~m 2 (140) 
which leads to: 

<™> 2 -?(S) 2=o (141) 

i.e. the principal symbol and the associated HJ equation do not take into account 
the mass of the particle. 

We see from these examples that the principal symbol does not capture the 
full physical information contained in the differential operator (with the associated 
PDE) we started from. To remedy this situation, we may define an associated, 
homogeneous, differential operator on a larger space by adding one more degree 
of freedom. To be specific, let us consider the case of a second-order differential 
operator. Locally, the operator will have the representation: 

D = a ^^- k +b ^ + e (U2) 

(aij,bj,e e Adding then one more variable, denoted as r, we obtain the 

extended differential operator on Q x M: 

„ ~ d d , d d d 2 , 

D - D = a ^^ k +b ^0-r +e 0^ (U3) 

which is now homogeneous of degree two. We will also restrict the space of functions 

on which D operates to functions of the form: 

f(x,r)=e T f(x) (144) 
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in such a way that: 

D (/) = e/D (/) (145) 

The polynomial function on T* (Q x IR) (with coordinates: (xj , r; p> , p T ) ) associated 
with the principal symbol of D will be: 

f 5 = a ]k pip k + b ]T Pp T + e (p T ) 2 (146) 

and we will look for solutions HJ equation: 

(ds)*(f 5 )=c (147) 

of the form: 

S(x,t) =t + S(x) (148) 
In this way, the HJ equation will become, when written in local coordinates: 

dS dS , dS 

and we have gotten rid in this way of the additional degree of freedom whose 
introduction was made necessary in order to be able to deal with tensorial objects 
(see the discussion in Sect. 3.1). 

Example 8. With the procedure outlined above, both the Schrodinger and the 
Klein- Gordon operators get replaced by: 

D ^ 5 ^' h ^Fr + ^ 2 - V ^& < 15 °» 

and: 

D KG - D KG = V 2 - IJ^ - (151) 

respectively, and the associated HJ equations (looking for solutions of the form 
(148)) will capture the full physics of the respective problems. 

The procedure outlined here is completely general, and we see that in this way 
the full information contained in our original differential operator will be captured 
by the principal symbol of the extended operator, and we can proceed now with the 
first-order PDE and Hamilton's equations just as before. As for transformations, 
we should remark that now we have to restrict to bundle automorphisms. 



3.3. Vector- Valued Differential Operators and the Hamilton-Jacobi Prob- 
lem. In several physical situations, the systems we may want to describe have also 
some inner structure, and therefore scalar differential operators are not enough. For 
instance, the quantum-mechanical description of a particle with spin requires that 
we replace Schrodinger's equation with Pauli's in the non-relativistic case and with 
Dirac' in the rclativistic case. More generally, this kind of situation occurs when 
we consider [8, 9] Yang-Mills fields and generalized Wong equations. 

Let us consider then the general aspects of this situation. We consider here two 
vector bundles: E\ — > Q and: E2 — ► Q- We denote by Sec(Ei) and Sec (£2) the 
spaces of sections: s : Q — > E\ and: r : Q — > E%- The operator of multiplication by 
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a function / £ J (Q) will be denoted here too by /. Following the algebraic setting 
for differential operators of Sect.3.1, we define: 

A differential operator of order at most k, acting from Sec(E\) to Sec{E2), is 
a linear map: 

D : Sec (E^ — ► Sec (E 2 ) (152) 

such that: 



D {k \h ,fi 



■j fk 



(153) 



V a ^ 



(154) 



Jk\ =0 V/ ,/i, 

In a chosen trivialization of the two bundles, their sections are vector-valued 
functions on Q, and the operator D will be described by a matrix whose entries are 
coordinate expressions of scalar differential operators. Thus, here too the relevant 
differential part can be written as: 

D (k) = 

\<j\<k 

a = (ii,i 2 , i n ) ; |er| = h + 12 + ■■■ + in 
d^/dx a = dW/dx\ 1 dx%...dx% 

but now, instead of the g a 's being pointwisc scalars, we have: 

g a e Hom(Ei,E 2 ) (155) 

(pointwise), and each matrix element of the g^s will be a smooth function on Q. 

Wc can again identify the principal symbol with a symmetric and totally con- 
travariant tensor, i.e.: 

d d d 



£>(*) 



dxi. 



dx,. 



(156) 



but now the tensor field must be understood as a multilinear function on one-forms 
with values in Horn (E\, E-i). Intrinsically: 

a p (D«) (dh,df 2 ,...,df k ) = [... ,.../ 2 ] ,/i] (157) 

and the symmetry follows again from Eq.(lll). By using the contraction with the 
fc-fold product of 9o with itself we obtain a polynomial function with matrix-valued 
coefficients. 

Example 9. Let; D = d, the exterior differential: d : A J (Q) — * A J+1 (Q). The 
value of its symbol on a differential one-form a is a homomorphism from A J (Q) to 
A i+1 (Q) given by [22]: 

a p (d) {a) : A J (Q) 3 uj i — > a A ui e A J+1 (Q) (158) 

The symbol of second-order scalar differential operators are symmetric con- 
travariant tensor fields of rank two, that is: 

d _d_ 

u 3 dx k ' 

When this tensor is not degenerate, it defines a (pseudo) Riemannian metric: 



(159) 



9D 



a 3 dxj ® dxk 



ith: 



a jk a km 



Si 



(160) 
(161) 
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Thus, D is elliptic if gD is Riemannian and hyperbolic if it is Lorentzian. 
In the scalar situation we have defined a Hamiltonian as: 

H(q,p)=a p (D)(0 o ®...®0o)= ,, fi9 v 
= a- ■ » 11 1 J 

— <-H 1 ...l k P ■■■P 

and in this way we have been able to define Poisson brackets, Hamilton- Jacobi-type 
equations and, finally, Hamilton's canonical equations. 

In the present context, the symbol of a matrix- valued differential operator D will 
be a matrix- valued polynomial on T*Q. For example, for a second-order operator 
we will get: 

a p (D){8 Q ®6 Q )=H(q,p) (163) 

where now H (q,p) will be a matrix. 

If our bundles were Hermitian bundles of the same dimension, they can be iden- 
tified and H becomes Hermitian: H = W. Then, H will have real eigenvalues, 
and each eigenvalue function can be used as a Hamiltonian function on T*Q [6]. 
At this point we can repeat whatever has been said in the case of scalar differential 
operators. 



4. Joint Hamilton-Jacobi Problems. As anticipated in the Introduction, one 
can pose a "joint Hamilton-Jacobi" problem also for two or more dynamical vari- 
ables. For reasons that will become apparent shortly, we will limit ourselves here 
to the "conventional" H J problem as discussed in Sect. 2. In order to discuss the 
"joint" problem, we will have to discuss first some preliminary notions related to the 
restriction of Poisson brackets (on a cotangent bundle) to Lagrangian submanifolds 
that are graphs of closed one-forms. 

In Sect. 2 we have shown that, given the cotangent bundle T*Q of a configuration 
manifold Q with the canonical symplectic form loq, with every one-form a £ X* (Q) 
we can associate the maps (52) and (55), i.e. (we omit here for brevity the suffix 
"a" that was employed in Sect. 2): 

p:Q^T*Qby:q^(q,a(q)) (164) 

and: 

ip =: ip o n : T*Q -> T*Q by : (q,p) h-> (q, a (q)) (165) 

and that the image of a: T[a] — ((q,a(q)) £ T*Q;q G Q) will be a Lagrangian 
submanifold of T*Q iff a is closed. Notice that, if (q,p) belongs already to T[a], 
the map (165) will leave it unaltered: 

^(r[a])=T[a] (166) 

in a fiber-preserving way. 

With every function / e 3 r (T*Q) we can associate the pull-back: 

(rf) (q,p) =/(«,«(?)) (167) 
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i.e. the extension of the restriction of / to T [a] with the property of being constant 
along the fibers. With this in mind, we want here to relate the restriction tp* {/, g} 
of the Poisson bracket of any two functions / and g to T [a] to the corresponding 
restrictions ip*f and ip* g of / and g. 

If 7r : T*Q — > Q is the canonical projection, Eq.(165) implies: n o ip = tt, and 
hence: 

Ttt o Tip = Ttt (168) 

for the corresponding tangent maps. 

Let now m = (q,p) £ T*Q and: 7 = ip(m). Then to every tangent vector 
X £ T m T*Q one can associate the vector: Tip (X) £ T 7 T*Q and Eq.(168) implies: 

T-ir(Tip(X)) -Ttt(X) = (169) 

Notice that the two arguments of Ttt in Eq.(169) are in general tangent vectors at 
different points of T*Q. It is only when m = 7 £ T [a] that we can factor out the 
tangent map Ttt and write: 

Ttt ((Tip (X)) -X)=0, Xe T-,T*Q, 7 £ T [a] (170) 

What Eq.(170) proves is that, if X £ T 7 T*Q, then (Tip (X)) - X is a vertical field. 

Let now /, g £ 5" (T*Q) and let Xf,X g be the corresponding Hamiltonian vector 
fields defined via: 

«x f ^o = -d/, «x s wo = -dg (171) 



while: 
Then: 



{f,g} = coo (X g ,X f ) (172) 



(173) 



(iP* {/, 5 }) (m) = {/, .9} (7) = ^0 (7) (* fl (7) , (7)) = 
lo q (7) (X ff (7) - Tip (X g (7)) , X/ (7) - I> (7))) + 
^0 (7) (* fl (7) , Ti> (X f (7))) + co (7) (Tip (X g (7)) , X f (7)) - 
loo (7) (T4>(X g (j)),T4>(X f (7))) 

Now, in the first term on the r.h.s. of this equation, too is evaluated on a pair of 
vertical fields (at 7), and hence this term vanishes. As to the last term, using the 
definition of the pull-back [41]: 

coo (7) {Tip (X g (7)) , Tip (X f (7))) = (ruoo) (7) (X g (7) , X f (7)) (174) 

But: ip*u>o = da and, as a is closed by assumption, this term vanishes as well and, 
using also Eq.(171), we are left with: 

{f,9}) M = (»T^(x / ( 7 ))ix 9 (7) - iT4>(x B (y))ix f (7)) ^0 (7) = 

(df) (7) (T^ (7))) ~ {dg) (7) (T^ (X/ (7))) - , 7 ,s 

(dr.f)(x g ( 1 ))-(dr 9 )(x g ( 1 ))= y 1 

(L Xg rf-L Xf rg) (7) 

The last term here is in turn the pull-back via ip of the function: Lx g tp* f — Lx f tp* g, 
and in this way we obtain the result [38] : 

r{f,g} = r{L Xg rf-L Xf 4'*g) (m 
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which expresses the relation between the pull-back (the restriction) of the Poisson 
bracket of any two functions and the restrictions of the functions themselves. All 
that has been proved here relics in a crucial way on T [a] being the graph in a 
cotangent bundle of a closed one-form on the base manifold, and cannot be extended 
straightforwardly [38] to more general contexts such as Lagrangian sumbanifolds in 
general symplectic manifolds without further qualifications (see however Ref. [38] for 
some possible generalizations). 

Having established the result (176), let's turn now to what we have called at 
the beginning of this Section the "joint HJ problem". To be specific, let: A,BG 
3 r (T*Q) be two dynamical variables. We look then for a closed one-form a such 
that (cfr. Eqs.(164) and (165)): 

ip* (A — a) = tp* (B — b) = 0; a, 6 e M (177) 

or: 

i'*A = a, ^*B = b (178) 
Eq.(176) tells us immediately that this implies: 

r {A, B} = r (Lx B ip* A - L XA rB) = (179) 
This [necessary) condition, i.e.: 

i/j*{A,B} = (180) 

can be interpreted [38] as a sort of Poisson theorem, i.e. a condition saying that 
if a solution exists for the joint HJ problem for A and B, then it must be also a 
solution for the HJ problem for the Poisson bracket {A, B} on its zero level-set, 
and also as a classical counterpart of the quantum condition [43] according to which 
two observables must commute in order to be simultaneously diagonalizablc. It is 
obvious that if the Poisson bracket {A, B} does not vanish anywhere, then there is 
no possible solution for the joint HJ problem for A and B. 

The obvious generalization to any set A\, Ak, k < n = dim (Q) will be, of 
course: 

r{A i ,A j } = 0;i,j = l,...,k (181) 

As the one-form a is closed by assumption, at least locally (and globally if the 
first de Rham cohomology group of Q vanishes): a = dW, W <G 5" (Q), and Eq.(180) 
will become a PDE for the unknown function W . For example, in local (Darboux) 
coordinates, cuo = dpi A dq l will imply the equation: 

~dA(q,p)dB{q,p) dB (q,p) dA{q,p) 
dq i dpi dq i dpi 

We will close this Section by discussing a simple example that can help clarifying 
the status of the condition (180) as a necessary but not sufficient condition. 

Example 10. It is well known [44] that the dynamics of the 2D harmonic 
oscillator onT*R 2 , with coordinates (q 1 , q 2 , Pi , P2) and equipped with the canonical 
symplectic form: loq = dpj A dq 3 can be described by in many alternative ways and, 
among others, by: 



Pi = mm = (182) 
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The "standard" Hamiltonian: 

H+ = Hi + H 2 : Hi = - 
and vector field: 



(ptf + fqf], i = 1,2 (183) 



r+ = r 1 + r 2 ; r t = Pi — - q l — j = 1,2 (184) 

oq l Opt 

or by: 

• The Hamiltonian and vector field: 

H_=H 1 - H 2 ; r_ = Ti — T 2 (185) 

or, eventually, by: 

H = p 1 p 2 + q 1 q 2 (186) 

and: 

r = p 2TT +pi-^ - q q — 9^— (187) 

ag 1 ag^ api op2 
which describes the dynamics of the 2D harmonic oscillator on the configura- 
tion space. 

Now, it is clear that: {H + ,H-} = 0, each Hamiltonian being a constant of 
the motion for the other, and the joint HJ problem for them is trivial, as complete 
integrals can be found simply by separation of variables, i.e. by setting: W(q , q 2 ) = 
Wi (q 1 ) + W 2 (q 2 ) , thereby splitting the HJ problem for H± into separate (and 
identical) one- dimensional HJ problems for H\ and H 2 . Considering instead the 
(alternative) Hamiltonian (186), we find again: {h+,h} = 0, but: 

{H_,H}=2(q 1 p 2 -q 2 p 1 ) (188) 

(i.e. the "canonical" angular momentum) and the PDE equation (182) becomes 
then: 

x dW ^ 2 dW 
dq 2 

whose only solutions 1 ® are of the form: 



W(q\q 2 )=W(r 2 ); r 2 = (q 1 ) 2 + (q 2 f (190) 

with W any (at least C l ) function. On the other hand, the HJ equation for the 
Hamiltonian H reads: 

dWdW , 2 . , 

ww +qq=E (191) 

and no function of the form (190) 2a can solve Eq.(191), nor it can for the HJ 
equation associated with the Hamiltonian (185). This clarifies the fact that the 
vanishing of the restriction of the Poisson bracket is only a necessary and by no 
means a sufficient condition for a solution of the joint HJ problem. 



19 Switching to plane polar coordinates (q 1 = r cos if>,q 2 = r sin 0) : q 1 d/dq 2 —q 2 d/dq 1 = d/d(p. 
20 Inserting (190) into (191) leads to the equation: q 1 q 2 (2(W') 2 + l) = E, which has obviously 
no solution. 
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5. Time-Dependent Hamilton-Jacobi Theory and the Hamilton-Jacobi 
Problem on Lie Groups. 

5.1. Time-Dependent Hamilton-Jacobi Theory. We will recall briefly here 
how the time-dependent H J theory as described, in local coordinates, by Eq.(l), 

16 rr( 9S \ OS „ 

H {*ir q > t ) + m= (192) 

can be recasted in a geometrical setting [38] just as its time-independent counter- 
part. 

Although the "minimal" extension of the carrier manifold [41] appropriate to the 
description of a time-dependent dynamics seems to be from T*Q to T*Q x R, where 
R stands for the time variable t, the resulting manifold, being odd-dimensional, 
is a contact [1] and not a symplcctic manifold. As anticipated in Sect. 1.5, it is 
more convenient [1, 38, 41] to extend the original configuration space Q (with local 
coordinates q' l ,i = 1, ...,n = dimQ) to: Q = Q x R, considering then time as an 
additional coordinate on the same footing as the g^'s, and hence to consider the 
(extended) tangent bundle: 

M =: T*Q = T*Q x T*R = M x T*R; M =T*Q (193) 

Coordinates 21 for T*R will be denoted as (t, h), with h the (energy) variable canon- 
ically conjugate to the time t, and a point fa G M will consist of a triple: fa = 
(m,t, h) with to G M (hence: m — (q,p) in local coordinates). Associated with the 
manifold M there will be various projection maps, for instance from T*Q x T*R 
to T*Q x R or to T*Q or to T*R and so on. One can further endow M with a 
canonical one-form 9 obtained by "adding" a contribution from T*R to 9$ = Pidq 1 
on M, i.e.: 

9 =: 0o + hdt (194) 

Hence, M will acquire the structure of a symplectic manifold with the canonical 
two- form 22 : 

uj = uj + dh A dt (195) 

The canonical equations of motion on M = T*Q can be "lifted" to M = T* (Q x R) 
by defining, guided by the structure of the time-dependent equation (192), the 
Hamiltonian H G ^(M) on the extended phase space M as: 

H = H + h (196) 

In this way: 

• The Hamiltonian vector fields Xh,X^ G X(M) associated respectively with 
H and H (i^^ujQ = —dH,i^_too = —dH) are easily found to be: 

OH d ~ ~ d 

Xh = Xh - -moE' x * = Xn + di (197) 



21 Global coordinates, in this case. 

22 To be precise, one should write on the r.h.s. of Eqs.(194) and (195) 7r*0o and tt'ujo, where 
it is the projection: n : T*Q X T*R — > T*Q, instead of 8o and a;o, but we will avoid, here and in 
what follows, explicit mention of the pull-back operations involved as long as this does not lead 
to ambiguities. 
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(X H = (dH/dp) d/dq - (dH/dq) d/dp), and hence: 
• Denoting with r the evolution parameter in M, the dynamics generated by 
Xfi is described by the set of ODE's in the parameter r: 

dq = 8H_ dp = _dH_ 
dr dp ' dr dq 

and: 

dt dh OH , s 

* = * = "7* (199) 
These equations imply: t — r = const, (and hence Eqs.(198) will "project 
down" to the canonical equations on T*Q). They imply also: dH/dr = 0, 
i.e.: H = H + h = const. 

Denoting with {., .} and {., .} the Poisson brackets on M = T*Q and M = 
T*(Q x R) respectively, we have: 

and, in particular, if: F G ^(T*^ x R) (i.e. it does not depend on h) and: 
G = H: 

{f : h}={F,H} + ^ (201) 

i.e: 

dF r ~i 

(202) 

and (possibly time-dependent) constants of the motion will be characterized 
by: 



{F.S} 



(203) 
in full analogy with the time-independent case. 

The dynamics on M leaves invariant (i.e. is tangent to) all the submanifolds 
where H takes a constant value, and the dynamics on T*Q x M can be recovered 
from any one of them. For this reason, one can concentrate on the zero-level set of 
H, i.e. on the invariant submanifold 23 : 

So = H^ 1 (0) = (to, t,h= -H (to, t)) C JVC; to € M, t G M (204) 

Being obtained by giving explicitly the (globally defined) variable h in terms of 
to and t, So is diffeomorphic to T*Q x R, and is actually a global section of the 
projection map: tt : T* (Q x R) -> T*Q x R. 

If we keep i fixed, we obtain a family {(Eo)t}tgK of constant-t sections of Eo- 
Each section is a submanifold of T* (Q x R) of codimension two (i.e. it is 2n- 
dimcnsional, while Eo is (2n + l)-dimensional) diffeomorphic to T*Q and is a global 
section of the projection map: tt' : T* (Q x 1) — > T*Q. This family provides 
then a (regular) foliation [41] of Eo- Moreover, due to the additional and nowhere 
vanishing d/dt which is present in the definition of X^, the (Eo)t's are transversal 
to the flow generated by X^, and the one-parameter group {$ r } rS R of canonical 



3 Being of codimension one, So will be a coisotropic [38] submanifold of OVt. 
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transformations of T* (Q x R) generated by will permute the (Eo)t's among 
themselves, i.e: 

$ T : (E ) t -> (£o)t+r; r G R; $ r o $ T , = $ r+T ,; $ = Jd^ (205) 

Remark 9. One should keep in mind that, while on the extended phase space 
T* (Q x 1), due to the fact that the dynamical vector field does not depend 
explicitly on the evolution parameter t, the dynamics is "autonomous" [1, 41], and 
gives rise therefore to a "bona fide" one-parameter group, this is not so on T*Q. 
There, due to the time-dependence, the dynamics will give rise in general [1, 38] to 
a two-parameter family {^t^} of canonical transformations of T*Q obeying the 
law of combination [38] : 

0*3*2 ° 0*2*1 = 0*3*1 ; 0*1*1 = Id T*Q (206) 

which will "collapse", of course, into a one-parameter group if the Hamiltonian 
happens to be time-independent, i.e., in such a case: </>t 2 ti = 4>t 2 -U .Explicitly (see 
also Eqs.(198) and (199)), the one-parameter group {$ T } will act as [38]: 

$ T : (m, t, h) -> (0 t+Tjt (m) , t + t, h + H (m, t) — H (0 t+T (m) , t + t)) (207) 

At this point, it becomes apparent how one can cast the time-dependent HJ 
problem as stated in Eq.(192) into a geometrical form on the extended phase space 
M, paralleling to some extent the discussion in Sect. 2. 2. With any S € ^{Q) 
{dS £ X*(Q)) we can associate the map: 

~ ~ ~ / dS dS\ ~ 

Ifs-Q-^M; (p 3 :QB (q,t) -> (q,—,t,—j CM (208) 

which is a global section w.r.t. the projection map: tt : M — > Q. Then: 

f =: <p 8 (Q) C M (209) 

will be the graph of the closed (actually exact) one-form dS on Q, and hence a 
({n + l)-dimensional) transversal Lagrangian submanifold in M. The requirement 
that S be a solution of the PDE (192) can be rephrased in geometrical terms by 
requiring that: 

tp* s H = <S=> f c S (210) 

Remark 10. Eq.(210) should make it clear why, at variance with the time-independent 
case, the "zero- level set" Eo plays a distinguished role in the geometrical approach 
to the time-dependent HJ problem. 



The geometrical meaning of the time-dependent HJ problem in the present for- 
mulation is therefore the following: 

We look for transversal Lagrangian submanifolds T in M that are contained in 
Eo and are graphs of exact one-forms on Q. T will have codimension n in Eo- A 
complete solution of the time- dependent HJ problem will be an n- dimensional foli- 
ation of the "zero-level set" Eo by such submanifolds. Usually, these submanifolds 
are defined by the "dispersion relations" of our PDE. 
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By reasoning as in Sect. 2. 2 we can conclude that the vector field will be 

tangent to the submanifold T. Hence, the dynamical flow will leave it invariant, 
i.e.: 

$ T (f)=fVT (211) 

All of the above picture is on the extended phase space M, and one should see 
now how it can be made to "descend" to the physical phase space M = T*Q, i.e. 
how it behaves under the projection: 

tt' : T* (Q x R) = M -> M =T*Q; tt' : (to, t, h) -> to (212) 

For every fixed t we may consider S as a function S t £ J (Q) with dS t £ X* (Q). 
Then, the map: 

<Ps,t ■ Q 9 Q (9, (dS t ) (q)) £ T*Q (213) 
giving a global section w.r.t. the projection: 7r : T*Q — > Q, will define the transver- 
sal Lagrangian submanifold T t = ip$ t (Q) m T*Q, the graph of the exact one-form 
dSt. 

Going back now to the situation in M = T* (Q x R), we can consider the family 
of intersections {T D (Eo)t}TeK of T with the constant-t sections of So- Explicitly: 

fn<£ ), = {(,,i|, f ,f)v t } (214) 

and, for every t, T n (£o)t w iU be an n-dimensional (hence isotropic [38] in M) 
submanifold contained in (So)*. Using then Eqs.(205) and (211) one sees that $ T 
permutes again these submanifolds among themselves, i.e.: 

$ r (f n(Eo) t ) =f n(s ) t+T (215) 

and, under the projection (212) we obtain: 

7r'(f n(So)t) =T t (216) 
Moreover, using Eq.(207), we see that the IYs evolve in time as: 

</H+r,t (r t ) = r t+r ; (,r£l (217) 
We obtain then the following picture: 

Any solution of the HJ problem in M = T* (Q X K), i.e. any single geometrical 
object r, gives rise in M = T*Q to a family {T t } of time-dependent, Lagrangian 
submanifolds in T*Q, each T t being the graph of an exact one- form dSt on Q, that 
evolves in time according to Eq.(217). 

One should keep in mind, however, that time-evolution need not be, and quite 
often is not, a harmless process, to the extent that solutions that are "well behaved" 
at a certain time may develop caustics and/or become ill-behaved at later times. 
To "cure" these and other possible pathologies, and also to fully implement sym- 
metries in the time-dependent context as well as in the time-independent one (see 
also the discussion at the end of Subsect.2.2), one might be forced to require that 
the relevant submanifolds be graphs of closed but not necessarily exact one-forms 
and/or to abandon the requirement of transversality. We will not discuss here these 
generalizations, but refer rather to the literature [38] for further details. 



HAMILTON-,! ACOBI THEORY 



35 



5.2. The Hamilton-Jacobi Problem on Lie Groups. In the previous Sub- 
section we have seen how, even for non-conservative systems, one can implement 
canonically, via the one-parameter group ($ T ) Te jj, the action of the Abclian group K 
of time translations by suitably enlarging the phase space from T*Q to T* (Q xl). 
It is interesting to generalize this approach [38]to the case in which the action of 
the Abelian group R is replaced by that of a more general Lie group G. 

5.2.1. Canonical Actions of Lie Groups on Symplectic Manifolds. Let M = T*Q 
with the canonical symplectic form 24 ojq and let g be the Lie algebra of a k- 
dimcnsional Lie group G, with a basis {e r } r=1 obeying the Lie bracket relations: 

[e r ,e s ]=C t rs e t (218) 

the C* s 's being the structure constants of the group. 

The group G will act canonically [1, 38, 41] on T*Q if there exists a realization 
{^g} s gG YIS " a f amnv °f diffeomorphisms of T*Q satisfying: 

<f>g O ff>g, = (f>gg, Vfl> , {/ G G <f> E = I (fa* Q (219) 

(with e the identity of the group) and: 

0*w o = w V.g G G (220) 

To the basis {e r } r _ 1 of the Lie algebra there will correspond a set {X r } r _ 1 , X r G 
X (T*Q) ,r = l,...,k of fundamental [41] vector fields generating the finite (canoni- 
cal) transformations (f> g and obeying the commutation relations: 

[X r ,X 8 ] = C l rs X t (221) 

As a consequence of Eq.(220) the X r 's will be locally Hamiltonian, i.e.: Lx r ^o = 
V?*. Wc shall assume them to be actually globally Hamiltonian, i.e. that there 
exist globally defined functions Hi,..., Hk s.t.: 

ix r tuo = -dH r , r=l, k (222) 

The Poisson brackets of the H's will be given as usual, by: {H r , H s } = ixJx r ^o, 
and it is a simple matter to show that Eqs.(221) and (222) imply: 

d[{H r ,H s } + C t rs H t ] = (223) 

i.e.: 

{H r , H s } = -C l rs H t + d rs ; d rs G R, d rs + d sr = (224) 
The Jacobi identity on the Poisson bracket (224) implies then: 

Cr S dut + cycl. perm, of (r, s,t) = (225) 

One can get rid of the additional constants if there exist other constants Ai, \k 
s.t.: 

C* s A t = d rs (226) 



24 Although we will concentrate here on T*Q, most of what will be said applies equally well 
[38] to a general symplectic manifold (JVC, u;). 
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in which case, re-defining: H r =: H r — X r one gets Poisson brackets without addi- 
tional additive constants, i.e.: 

{H r ,H s } = C^Ht (227) 



Remark 11. The condition (225) is the statement that the map: d : g x g — > E 
acting via: d(e r ,e s ) =: d rs (with g acting trivially on R) must be a two-cocycle 
[17, 37, 39, 40] in g, while Eq.(226) requires it to be also a coboundary. This will 
be granted if the second cohomology group H 2 (g,R) of g with values in M [17] 
vanishes, i.e.: H 2 (g, M) = 0. This is true for many of the groups that are relevant 
for Classical Dynamics [51] as well as for Quantum Mechanics and Field Theory such 
as the Euclidean, the Lorentz and the Poincare groups. It fails however to be true 
for the Galilei group, which has a non-trivial two-cocycle (i.e. a cocycle which is 
not a coboundary) connected with the mass [51] which is an invariant for the Galilei 
group. 

5.2.2. The Canonical Action of a Lie group G on T*G. As is well known [41], a Lie 
group G can act on itself, among other ways, via left or right translations. We will 
concentrate on the former, i.e. on the diffeomorphisms: 

Lf : g> - 4°V =: 9sf\ 4? L fJ = L fL'> V<7, </,<?!, <?2 G G (228) 

and will discuss briefly how this action can be lifted to a canonical action of G on 
the cotangent bundle T*G. 

Left translations are generated by a set : x[°\ of independent and nowhere 

vanishing right-invariant [41] vector fields satisfying: 

with associated dual (right-invariant) one-forms 9^ r ,r = 1, k such that: 

i x i 0) 6 mr = 5l (230) 
and satisfying the Maurer-Cartan structure equations [41]: 

de (o)r + I(7j t 6K°) S a 6»( 0)t = (231) 

These right-invariant one-forms will provide a basis 9^°' r (g) for the cotangent space 
T*G at any point jeG. 

A basis for the dual g* of the Lie algebra g will be denoted as {e r }* =1 and will 
be specified by: 

(e r \e s } = 6; (232) 

Any element h G T*G can be written as: h = h r 9^ r (<?), and we have the associa- 
tion: 

T g *G 3h = h r 6 {a)r (g) <^h = h r e r G 0* (233) 
This implies [41] that T*G will be diffeomorphic to G x g*: T*G ~ G x g*. 



(0) 



(229) 
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Using the canonical projection: 7r : T*G — ► G we can pull back the 9^ r, s to 
obtain: 9 r =: Tr*9 (Q ^ r e X* (T*G) ,r=l, k, and we can define then the canonical 
one-form Go on T*G as: 

6 = h r 9 r (234) 

Hence: Slo = d<do will be the canonical symplectic form on T*G. 

The canonical lifts X r 0) * from G to T*G of the 's will be defined [41] by 
the conditions: 

L x ( 0) ,e = 0, Tir (.X^*) = r = 1, .., k (235) 

i.e. Xr must be a vector field that leaves the canonical one-form unchanged and 
that projects down under the tangent map T7r to X^ for all r's. Xr is of course 
Hamiltonian as the first of Eqs.(235) is equivalent to: i„(oj»fio = — e£(i„(o)»Oo). The 
X^^'s turn out also [38] to be dual to the 6 r, s. i.e.: 

i x m*e s = 5 s r ^i x m,Q = h r (236) 

and hence: 

i x p).£lo = —dh r (237) 

The Xr ^*'s close on the same Lie algebra as the X^'s [51], i.e.: [Xr " 1 *, -^i '*] = 
C* s X t (0) *. The same happens with the Poisson brackets of the /i r 's, but in this case 
(cfr. Eq.(224)) without [41, 51] additional constants i.e.: 

{h r ,h s } = -C t rs h t (238) 

with no need to require here the second cohomology group iJ 2 (g,R) to vanish. 

It is the lifted vector fields Xr that generate the canonical action of G on T*G. 
Finite transformations will be denoted as T*L^ and they will act as [38, 51]: 

T*Lf) : (g',h') — (gg' , V (g)h') (239) 

where D (g) is the coadjoint [41] representation of G on g*. 

5.2.3. The Hamilton- Jacobi Problem on G. We shall extend now our formulation of 
the TDHJ problem by replacing K with G. We can combine the canonical actions 
of G on T*Q and on T*G ~ Gxg* to get a canonical action on: JVC =: T* (Q x G). 
M will be a symplectic manfold with the canonical two-form: 

uj = LO + Q (240) 

Points in M will be denoted as: fh = (m,g,h),m £ T*Q,g £ G and h € g* . The 
required canonical action will be denoted as <& g ,g € G, and will be given by: 

$ g : (m, </,/?) = (^(mW,D(ff)fc') (241) 
It will be generated by the vector fields: 

X r = X r + X^*,r = 1, k (242) 
which will obey the commutation relations: 

X r ,x}=C t rs X t (243) 
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and will be Hamiltonian w.r.t. the symplectic form u> with Hamiltonians: 

H r = H r + h r (244) 

One can prove, just as in Subsect.5.2.1, that the Poisson brackets of the H r 's are 
given by: 

{Hr,H a }=-Ci. B 5t + d ra (245) 

and they will "inherit" the same two-cocycle d that we found to be present there. 
We can define now the momentum map [41] \x as: 

fx : M — > q* : \x (m) = H r (m) e r = (H r (m) + h r ) e r (246) 

and define next the analog S of the submanifold Eo of Subsect.5.1 will be defined 
as the zero- level set of the momentum map, i.e.: 

E = ir 1 (0) = {(m,g,h\ \m e T*Q,ge G,h r = -H r (m) »F r = oJcM 

(247) 

which is a submanifold of dimension (2n + k) diffcomorphic to T*Q x G. 
Now, as: 

[H r , H S } = L ^ H r = -C l rs H t + d rs (248) 

we see that S will be invariant under the canonical action $ g if and only if d = 0. 
This is one of the most compelling reasons for requiring d to be a coboundary, i.e. 
for requiring that H 2 = 0. 

In analogy with the constant-t sections introduced in Subsect.5.1 we can intro- 
duce "constant-g sections" S 9 : 

E s = Um,g,h)\m G T*Q, g fixed, H r = oj (249) 

which are of dimension In and diffcomorphic to T*Q. The canonical action $ g 
permutes these sections among themselves, i.e.: 

% (? 9 >) = ^ 99 ' (250) 

One can pose now a geometrical H J problem for the Lie group G as follows: 
Find a (maximal) Lagrangian (hence oj dimension (n + k)) submanifold T C 
M = T* (Q x G), transversal w.r.t. the projection: tt : M — > M and such that: 

i~fj, = 0<* i~H r = 0, r = 1, .., fc => f c S (251) 

and sucht that, for every g £ G: 

r 9 = tt (f n e 9 ) =: r g (252) 

is a g-dependent Lagrangian submanifold in M = T*Q. 
This implies, of course: 

$ 9 (f ) = f (253) 

as well as, at the level of M: 

9 (IV) = IV (254) 
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Remark 12. Requiring the simultaneous vanishing of all the H r 's on V is remi- 
niscent of the joint HJ problems that were discussed in Sect. 4. The vanishing of 
the two-cocycle d is then another necessary condition for the HJ problem on Lie 
groups to admit of a solution. 

Finally, if L happens to be the graph of an exact one-form dS £ X (Q x G) 
determined by 5 € J(Q x G), we can express everything in somewhat more familiar 
terms. 

In the notations of Refs.[38, 51], let: a = (a 1 , a fe ) be local coordinates for G 
in a neighborhood of the identity, normalized to a r = 0, r = 1, k, at the identity. 
Then, if a, (3, 7 are coordinates for a, b, c G G, the composition law: ab = c is 
expressed in local coordinates as: 

7 r = f r (a,p), r=l,...,k (255) 

with prcassigncd functions 25 /''. Introducing the (invcrtible [51]) matrix: 

r,(a) = K (a)|| ;t£ (a) = ^^'^ l^o (256) 

and its inverse: £ (a) = ?y (a) 1 (yf s (0) = (0) = <5£), right-invariant vector fields 
are given [38]by: 

X(°> (<*) = -# (a) (257) 
and the associated right-invariant one-forms will be given by: 

0(O)r = _£r da s ( 258 ) 

The canonical one-form 0o of Eq.(234) can then be rewritten as: 

So = -h s &da r =: ir r da r (259) 

and will define the " momenta" : 

7r. r = —£*h a O h r = -if r ir s (260) 

Then, with: S — S (q; a) and the usual replacement: ir — > dS/da, the conditions: 
H r = H r + h r = become: 

/ dS(q; a) \ dS (q;a) 

Hr [ q > ~^q— ) ~ * ^ = ° (261) 

and this is now the (conventional) system of PDE's associated with a canonical 
realization of the Lie group G on T*Q. 

For further generalizations we refer once again to Ref . [38] . 



5 The conditions to which the / r 's have to obey are discussed at length in Ref. [51]. 
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6. Concluding Remarks. Following Dirac's prescription [21] according to which 
Classical Mechanics must be a suitable limit of Quantum Mechanics, we have con- 
sidered in this paper the Hamilton- Jacobi theory as emerging from Quantum Me- 
chanics when we consider an approximation in which Planck's constant is treated 
as a parameter. 

^From this point of view, while we keep on with the originating idea of the 
theory of considering Hamiltonian Optics as a suitable limit of Wave Optics, we 
have taken advantage of the fact that Wave Mechanics concerns also with particles 
with internal structure. This has suggested that we deal not only with "scalar 
differential operators" like Schrodinger's or Klein-Gordon's, but also with "matrix- 
valued differential operators" as those appearing in the Pauli as well as in the Dirac 
equations. The net result is an extension of the usual Hamilton- Jacobi formalism to 
a formalism where the Hamiltonian " scalar function" is replaced by a matrix- valued 
Hamiltonian. 

In the same spirit, following what happens in relativistic field theories, we have 
replaced the one-parameter group of time evolution with a Lie group, e.g. the 
Poincare group. 

Wc have not considered the Hamilton- Jacobi theory for field theories proper but, 
following the ideas outlined in this paper, it should not be very difficult to foresee 
how to proceed. 

The joint Hamilton-Jacobi problem turns out also to be very useful to deal with 
holonomic and non-holonomic constraints within the Hamilton-Jacobi formalism 
[15]. In this connection wc should also mention some recent relevant contributions 
[28, 34] to the same subject. 

The geometrization of the Hamilton-Jacobi problem presented here has many ad- 
vantages over more conventional presentations. For example, as mentioned already 
in Sect. 2. 2, posing a "Hamilton-Jacobi problem" as the search of foliations of the 
energy surfaces by Lagrangian but not necessarily transversal submanifolds opens 
the possibility of implementing the full set of canonical symmetries as symmetries 
of the Hamilton-Jacobi problem as well. We have also shown elsewhere [38] how to 
deal, in the same fully geometric spirit, with transformations and symmetries for 
partial differential equations of the Monge- Ampere type. Using the present gener- 
alization to differential operators acting on sections of vector bundles it should be 
possible to incorporate into the formalism more general PDE's than those of the 
above " Monge- Ampere" type. 

Using our matrix-valued Hamiltonians it will be possible to deal with equations 
of the Wong type [8], i.e. equations describing particles interacting with Yang-Mills 
fields. Also, in this approach, treating, say, electrons moving in some monopolc-likc 
magnetic field, algebroids arise in quite a natural way. We shall postpone more 
details on these aspects to a forthcoming paper. 
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